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ABSTRACT 

It is proved that for a system of spins <7j = ±1 having an interaction energy 
— KijO~iO~j with all the K\a strictly positive, one can construct a dual formulation 
by associating a dual spin S{jk = ±1 to each triplet of distinct sites i,j and 
k. The dual interaction energy reads —^2uj\DijYik^ijSijk with tanh(Kij) = 
exp(—2Dij), and it is invariant under local symmetries. We discuss the gauge- 
fixing procedure, identities relating averages of order and disorder variables and 
representations of various quantities as integrals over Grassmann variables. The 
relevance of these results for Polyakov's approach of the 3D Ising model is briefly 
discussed. 



1. Introduction, Basic Ideas and Main Results 



The order-disorder duality ' ' has been a useful tool to approach statistical 
or field theory models. In particular, the self-duality of the two-dimensional Ising 
model with nearest neighbor interactions on a square lattice provides an easy way to 
obtain the critical temperature and suggests a convenient way to write the partition 
function and the specific heat. The model dual to the three-dimensional nearest 
neighbor Ising model is a gauge theory' 131 and interesting loop equations, combining 
the order and disorder variables, have been obtained by Polyakov! 111 These equations 
suggest that the model near criticality can be described in terms of a string model. 

[31 

Despite interesting results, several aspects of this idea await clarification. 

Recently, we have proposed to study related problems in the case of the hierar- 
chical model. For this model, the renormalization group procedure can be reduced 

[2] 

to a local recursion relation which has been studied in detail and rigorous results 
concerning its fixed points and the spectrum of its tangent map have been ob- 
tained. This recursion relation is closely related to the the "approximate recursion 
formula" used by Wilson 1141 to estimate (rather well) the critical exponents of the 
three-dimensional Ising model. This motivated us to try to express the renormal- 
ization group transformations in various reformulations of the hierarchical model. 
In particular, we constructed the dual formulation and the several representations 
of the partition function in terms of Grassmann variables. Due to the non-locality 
of the model - the price to pay for the locality of the renormalization group trans- 
formation - appropriate methods are necessary to enumerate the terms appearing 
in the high (or low) temperature expansion. The methods we found apply to a 
larger class of Ising models than the one considered initially and are reported here. 

In this paper, we consider the class of systems of spins U{ = ±1, where % runs 
from 1 to N, with an interaction energy 



and all the iV(iV — l)/2 couplings Kij strictly positive but otherwise unspeci- 




(1) 



i<j 
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fied. This class of models includes for instance, the hierarchical model and one- 
dimensional Ising models with \i — j\~ a interactions. 



As usual, we define 



and 



< a. 



n ■ 



E ( 2 ) 

{<n=±i} 



.a i2n >=Z- 1 £ a h <J i2n e-P H . (3) 

{<T<=±1} 



In the following, a sequence of indices is understood as strictly ordered unless 
specified. For instance, the notation Kij implies % < j. We use parenthesis when 
the ordering needs to be performed. For instance, Kuj\ is equal to K^ if % < j 
and to Kji if j < i. This ordering is a matter of notational convenience and has 
no intrinsic significance in the problem. We also define the dual couplings by the 
relation 

tanhKij = e" 2Aj . (4) 

It is clear that the strict positivity of Kij implies the strict positivity of Dij. 

One of the main results obtained here is a dual formulation of the class of 
models defined above. The disorder variables are attached to each triple of distinct 
sites ijk and denoted Se^y The interaction energy can be constructed out of 
products of S(yfc) with i and j fixed and k running over all its values except for i 
and j. More precisely: 

Theorem 1 . The class of Ising models defined in Eqs. (1-4) admits the following 
equivalent (dual) representation 



___ _ +EP«( n S (ijk) -i)) 

Z = ([[coshKij) 2 N -y 3 j ^ e ^ k ^ (5) 

i<j {S I3fc =±l} 



3 



and if all the pairs ii^,...., %2n-\in ore distinct 




-Di 



«2n-l°"«2n 



*2n-l*2n 



> 



(6) 



In addition, for any four distinct sites i,j,k and I , the dual interaction energy is 



This theorem is proven in three steps in the next section. The first is a standard 
high-temperature expansion. The second establishes a connection between the 
terms of this expansion and some elements of a vector space. The third uses the 
fact that the homology of a certain boundary operator over this vector space is 
trivial to rewrite Eqs. (2) and (3) under the form given in Eqs. (5) and (6). This 
proof makes the local invariance mentioned at the end of the above theorem rather 
transparent. It is possible to eliminate this degeneracy by imposing an appropriate 
gauge-fixing condition. An example is given in section 3. 

Theorem 1 gives a procedure to express average values of the order variables in 
terms of average values of the disorder variables. This procedure can be inverted 
using the standard character expansion for Ising spins. Simple identities obtained 
this way are shown in section 4. In the case of short range Ising models, identi- 
ties relating appropriately chosen products of order and disorder variables can be 
interpreted as a discrete version of some kind of Dirac (or Ramond) equation 111 ' 31 . 
These equations are closely related to the Schwinger-Dyson equations associated 
with a representation of the partition function as an integral 1121 over (anticommut- 
ing) Grassmann variables.' 11 In section 5, we rewrite the partition function of Eq. 
(2) as an integral over Grassmann variables (Theorem 2). This formulation is very 
compact and the Schwinger-Dyson equations follow straightforwardly. Note that 
the integral representation is not minimal, it involves twice as many variables as we 
normally need. The elimination of these auxilliary variables is discussed at the end 
of the section 6 where Theorem 2 is proven. In the conclusions, we briefly explain 
how these results can be used in a treatment of the hierarchical model inspired by 
the work of Polyakov, Dotsenko and other authors. 



invariant under the simultaneous changes of sign of S^jfy, S^jfy, S^i) and S^iy 
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2. Proof of Theorem 1 



The first step in the construction of the dual formulation is a standard high- 
temperature expansion of the partition function and the correlation functions. We 
shall recall it briefly in order to introduce useful notation. We use the character 
expansion exp(Kij<Ji<jj) = cosh(Kij) + sinh(Kij)ai<7j and we introduce at each 
link ij a new variable nij which takes the values or 1. These can be used to 
rewrite Eqs. (2) and (3) through the identity 

ri( i +^)= e n4? w 

i<3 {n tJ =0,l} i<j 

The non-zero contributions come from the terms where all the o~i appear an even 
number of times. The high temperature expansion can thus be written as a sum 
over paths were the links ij appear at most once and the number of "visits" at the 
site % 

i-i N 

N i = E n ^ + E n y = E n (y) ( 8 ) 

is odd if i is one of sites of the spins appearing in the correlation functions and 
even otherwise. In the following we shall use a = 6(2) as a short notation for 
a = b modulo 2. The above arguments can be summarized by the first lemma. 

Lemma 1 For the class of Ising models defined in Eqs. (1-4) 

E °ii v^e-f> H = ([[coahKij) 2 N E e~ 2 5 A ^ (9) 

{cr; = ±l} i<j I «ij=0,l | 

1 Ni=mi(2) J 

with mi = \ if i G ,^2n} and otherwise. 

The second step consists in establishing a correspondence between the terms 
of the above expansion and the elements of a vector space over the field of integers 
modulo 2. Our goal is to express the condition Ni = mj(2) appearing in Lemma 
1 in terms of a boundary operator having a trivial homology. For this purpose we 
introduce the following definitions. 
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Definitions. Let V be a 2 N dimensional vector space over the field of integers 

modulo 2. To any sequence of integers 1 < %\ < i 2 < < i n < N, we 

associate an element of a basis of V denoted \i1i2 i n >■ The elements of this 

basis corresponding to sequences with exactly n indices span a subspace denoted 
V n . We define the boundary operator d n : V n — > V n -\ for 1 < n < N such that 

dn\hl2 in >= |u«2 in > +|*1*2 in > + + |*1«2 in > (10) 

where the hat designates a skipped index. We also define D n : V n — > V^+i for 
< n < — 1 such that 

n ,. • • ^ / I 1 ^1^2 in > if ii > 1 

-Dn|^2 l n >= < (11) 

I otherwise . 

The choice of the first index is as inessential as the ordering itself. 

We can now establish a one-to-one correspondence between any choice of con- 
figuration {riij} appearing in the high temperature expansion (9) and the vector 
of V 2 

\{mj} > =^2mj\ij > . (12. a) 

i<j 

Similarly, we can denote the elements of V\ 

\{m % } > = ^mi\i > . (12.6) 

i 

From (10), c^Kj >= \i > + \j >■ This relates in a clear way the link ij and the 
two sites "visited" % and j. Collecting all the terms and using Eq.(7) we obtain the 
following. 

Lemma 2. The condition Ni = mj(2) can be rewritten as 

d 2 \{nij} > = \{ mi } > . (13) 

In order to find the general solution of Eq.(13), it is crucial to note that the 
homology of d n is trivial. In order to establish this result, we first prove the 
following. 
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Proposition 1. Ifl<n<N — l, then 



a) d n d n+1 = 

b) D n -id n + <9 n+ iA 

c) D n+1 D n = . 



1 



(14) 



This follows in a elementary way from the definitions (Eqs. (10)-(11)). First, <9 n +l 
deletes one index in all the possible ways and d n d n+ \ deletes all possible pairs of 
indices twice. Since V is a vector space over the integer modulo 2, this proves a). 
Similarly, we obtain b) and c) as a mere application of the definitions by discussing 
separately the cases where the first index is 1 or is not 1. 

Corollary 1. // 1 < n < N — 1, then Im d n+ i = Ker d n . 

Part a) of the Proposition 1 implies that Im d n+ i C Ker d n . If \v n > G Ker d n 
then part b) implies that \v n >= d n +iD n \v n > which proves that Im d n +i D 
Ker d„. 



which amounts to the second corollary. 

We are now ready to rewrite restricted sums like those appearing in Eq (9) as 
unrestricted sums. 

Lemma 3. /// is a function V2 — > R and \ {mi} > e V\ then 



Corollary 2. Dim Im d n+ \ = ( N n 1 )- 

Using that Dim V n = (^) = Dim Im d n + Dim Ker d n and Corollary 1, we 
find Dim Im d n +i = (^) — Dim Im d n . Repeating n times we obtain, 




(15) 



E /(IK}>) = 2-(V) f(\Wj}> +ds\vs>) (16) 



{ 
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where \{n°j} > is any solution of 

d 2 \{n° j }>=\{m i }> (17) 

Two arbitrary solutions of Eq.(13) differ by an element of Ker 82 which is equal 
to Im 83 by Corollary 1. Consequently, we can write any solution as the sum of as 
particular solution > plus ds\v^ > for some \v% >G V3. This element of V3 

is determined up to an element of Ker d% = Im 84, a subspace having dimension 
( N ^ 1 ) by Corollary 2. In other words, if we consider now / to be a function over 
V3, it is invariant under the "gauge transformation" 

\v 3 >-> \v 3 > + d 4 \v 4 > (18) 

with |t>4 >G V4. When taking the unrestricted sum over V3, we need to divide 

(N-l\ 

by the "gauge multiplicity" which from Corollary 2 is 2v 3 ) , This concludes the 
proof of Lemma 3. 

We are now in position to prove Theorem 1. Using Lemma 3, we can replace the 
restricted sum over V2 (i.e., over the riy) appearing in Lemma 1, by an unrestricted 
sum over V3. More precisely, for any triples of distinct ordered sites i,j and k, we 
introduce a new variable taking the values or 1. A given choice of 
defines an element of V3 in an obvious generalization of Eqs. (12). We need to 

replace n-ij by (n°- + Yl <%jfc)) m °dulo 2 considered as a real number taking the 

k;k^i,j 

values or 1. This can be conveniently done with the substitution 

2 nij -> 1 + (-1) . (19) 

If the mi are like in Lemma 1, we can choose for instance, n? = = n? s = 

1 the others being zero, as in the particular solution appearing in Lemma 3, 
provided that all the pairs 1112,...., %2n-\in are distinct. The local invariance of 
the new formulation has been made clear in Eq.(18). Defining for convenience 
Sijk = {— l) diJk , we obtain Theorem 1 as stated in the introduction. 
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It is instructive to check that the the original model can be recovered if the 
duality transformation is applied to the dual model. Starting from Eq.(5) one 
can proceed as above and prove that the cohomology of a coboundary operator 
(algebraically dual to the boundary operator considered above) is trivial. The 
original model is then recovered easily. Since no new concepts or results are involved 
in this proof, it has been left as an exercise. 

3. An Example of Gauge-Fixing in the Dual Formulation 

The proof of Lemma 3 has made clear the fact that the sum over V 3 can indeed 
be rewritten as a sum over V 3 /Im 84. In this section, we show that the gauge 
condition 

D 3 \v 3 > = (20) 

is sufficient to fix a unique representative for each class of equivalence of V 3 /Im 84. 
This choice is motivated by the fact that we have studied the properties of D 3 
in the previous section. It is not a "covariant" choice since the first site plays a 
special role. When studying particular models, existing symmetries may suggest 
more appropriate choices. 

In order to show that the condition (20) specifies one and only one element of 
V 3 among those differing by a gauge transformation (18), we introduce a projection 
operator P defined as 

P = 1 + d4D 3 (21) 
From Proposition 1, it is immediate that P 2 = P and more importantly that 

D 3 P = Pd 4 = . (22) 

Given any element \v 3 > of V 3 , this implies that P\v 3 > satisfies the condition (20) 
and that if two elements differ by a transformation (18) their images under P are 
identical. This concludes the proof of the statement made above. 
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In terms of the coefficients in the \ijk > basis, the condition (20) amounts 

to 

d ijk = 0ifl<i<j<k (23) 

the remaining ones (dijk) being unspecified by the condition. It is easy to check 
that a similar result may be obtained starting with the high-temperature expansion 
(9) and eliminating the ny using the N — 1 independent restrictions on the iVj. 

N 

(It is clear that £ N { = 0(2).) 

i=i 

In general, a gauge condition of the form G\v% >= specifies uniquely the 
representatives and can be enforced by a projection operator Pq = 1 + d^H^, 
where H$ is an operator from V3 to V4 depending on G, provided that we can 
satisfy the consistency conditions G = Gd^H^ and 84 = d^H^d^. 

4. Simple Applications of Theorem 1 

In this section, we work out simple applications of Theorem 1. We explain how 
to introduce order variables in the dual formulation and disorder variables in the 
original formulation. We define the average value of a function of the dual spins 
Sijk as the r.h.s of Eq.(5) with this function inserted in the sum over the dual 
spins, divided by Z. We shall use the same brackets as in Eq.(3), their content 
indicating unambiguously which kind of average is considered. 

We know from Theorem 1 that an insertion of a^aj in an average value amounts 
to a change Dij — > —D\j in the dual formulation. This change can be implemented 
by an insertion of exp(—2D{j FJ S^j^) provided it occurs only once. The sim- 

plest application of this procedure is 

-2Dij Yl S (ijk) 

< oioj > = < e > ^ 

> . 

The positivity of the average follows from the positivity of the exponential and 
yields the first Griffiths inequality. 141 The generalization to an arbitrary number of 
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distinct pairs is straightforward. An interesting example is 



< (cri(Tj)(aja k )(a k ai) > = < e 




= 1 



(25) 



The second equality can be seen either from the fact that the square of the sigmas 
is 1 or from the invariance of Z, as written in Eq.(5), under the change of sign of 



The above examples show how to express products of an even number of order 
variables in the dual formulation. The procedure can easily be inverted in order 
to express gauge-invariant products of the disorder variables within the original 
formulation. The practical implementation is that exp(—2Kij<Ji<jj) corresponds 
to | [ S^-fc) in the dual formulation. This follows from a character expan- 

k;k=£i,j 

sion of the exponential, the expression of OiOj in the dual formulation discussed 
above, another character expansion and the identities sinh2Dij sinh2Kij = 1 and 
cosh2Dij sinh2K{j = cosh2Kij. Obviously, this implies 




< 



II s m > =<e 




> . 



(26) 



The analog of Eq.(25) reads 



< 



n ( n %^) > = < e 



-2 J2 KijcncTj 



> 



k;k^i,j 



(27) 



= 1 



Note that each SV^-m appears twice in the above products. Using similar methods, 



we obtain easily 




m;m^i,j n;n^=k,m m;m^i,j n;n^k,m 



(28) 



sinhKij sinhKki(< ^aja^ai > — < <Ji<jj >< o^oi >) . 
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The positivity of the r.h.s is the second Griffiths inequality! 41 

The above procedures can be used to write products of order and disorder vari- 
ables in both formulations. In the case of short range Ising models, it is possible 
to write suggestive identities relating average values of suitably chosen products of 
this type. A better insight concerning these equations may be obtained by reex- 
pressing the partition function as an integral over Grassmann variables. Note that 
this representation seems also well suited for the renormalization group method' 101 

5. A Representation of Z as an Integral over Grassmann Variables 

In this section, we describe a representation of the partition function as an 
integral over Grassmann variables (Theorem 2). We use this integral representation 
to express average values of the type considered in the previous section and to 
obtain Schwinger-Dyson equations. The proof of Theorem 2 is given in the next 
section. 

The basic properties of a set {^ a } of Grassmann variables can be summarized 
as follows. 

^1 =0 
m a ^ b = - m b ^ a 
d^ a m b = - m b dV b (2g) 

d^> a =0 



/ 



/ 



dm a m a =1 . 



A detailed presentation on this subject can be found in Berezin's book.^This tech- 

[12] 

nique has been used previously for short-range Ising models. In the following, 
Grassmann variables are denoted ip or \ with various set of indices. 

For notational convenience, we define 

cosh(x) if N is odd ; 
z{x) = { (30) 
sinh{x) if iV is even 
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/(*) = —ln(z(x)) (31) 



[diPdx] = Hdxid^dxid^ (32) 

i<j 

The following theorem will be proven inthe next section 

Theorem 2 . The partition function of the Ising models defined in Eqs. (1-4) 
admits the integral representation 

Z=C[[ ^shK tj )2 N J [dijdx] e ^ J] z{ £ (33) 

i<j ' i=l 

As in the previous section, we define the average of a function of the Grassmann 
variables by inserting this function in the integral representation of Z given in 
Theorem 2 and dividing by Z. Again we use the same brackets, their content 
raising the ambiguity. It is easy to express the average value of products of order 
or disorder variables considered in the previous section in terms of average values 
of functions of ipjipj- In particular, comparing the derivatives of Z with respect to 
Kij in Eq.(2) and Theorem 2, we obtain for i < j 

< OiOj > = thKij + (coshKij)" 2 < tpjipj > . (34) 

On the other hand, by changing the sign of in these two expressions of Z, it 
follows that 

< e - 2Kij<Ti<T i > = 1 - thK^ < ^) > . (35) 

These operations can be repeated in an obvious way, in order to obtain more 
involved products, provided that distinct pairs of indices are used. 

It is also possible to use the well-known fact that a constant shift does not 
change an integral over a Grassmann variables to obtain useful identities. This type 
of identities are often called equations of motions or Schwinger- Dyson equations. 
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Let A be an arbitrary function of the Grassmann variables and \& either ipj or 
x\ with i < j. Using the invariance of < A > under a shift in ^, we obtain the 
identity 

< {^KthKiiWiti -Xfrj] +/(£ #xJ)^[£ #X?]) A + >= ° (36) 

The function /(x) has been defined in Eq.(31). In the case N even, coth(x) must 
be multiplied the factor sinh(x) appearing in the integral before being Taylor 
expanded. We have used the usual convention that the derivative with respect to 
^ anticommutes with the other Grassmann variables. As a simple example, for 
A = ipj and \I/ = ipj the identity reads 

< > = (thK tJ )-\l+ <fC£ $x\)xWi >) • (37) 
Similarly, if \1/ takes the values ipj or x) with i < j, we obtain 

< - X{X)\ +fC£ ^4)^E A A >= o (38) 



6. Proof of Theorem 2 and Comments 

In this section, we prove Theorem 2, we discuss the extension of the method 
to the dual case and we discuss the elimination of the "auxiliary" variables x- 

The proof of Theorem 2 is based on the following identity 

1 + thKij(Ti<T 



I 



ij u i u j 



(39) 



E ("O n *Wx}) (1 - ntf) E (^r^jx)) 



j h(l-mij) 



71^=0,1 mij=0,l 



After using this identity for each link, we can rearrange inside the integral all the 
terms having a common subindex i, i.e., a common a{. Summing over the <7j, we 
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obtain a "local" condition namely that 

i-l N 

E m ^ + E n ^ = ( 2 ) ( 4 °) 

By construction, the and rnj do not appear at other sites and their sum can 
be performed locally. When N is odd (even), we obtain the sum over all products 
of an even (odd) number of distinct ipfxj- Using the fact that the square of a 
Grassmann number is zero, we can express this sum as the cosh (sink) of the sum 
over j of the ipjxf- This concludes the proof of Theorem 2. 

The dual version of this construction is easy to obtain. To each link ij we 
assign the Grassmann variables ipfj and x% with k running over all possible values 
but i and j. The non-local part of the integrand can be chosen as 

Y.thDij n v£ te n x% 

g i<j k\k^i,j g i<j k\k^i,j (41) 

(iV-2)(jV-3) 

where the products are ordered in the natural order and 7 = (— 1) 2 . Pro- 
ceding as in the above proof, we obtain the local part of the integrand 

II sinh ^ijXij + rtkXi k + i>Ux%) ■ (42) 

i<j<k 

The simplicity of the proof of Theorem 2 comes from the fact that we were 
able to rearrange the tpj \\ into local products. This was the reason to introduce 
the auxilliary variables x{ which "screen" the anticommuting nature of ipj. Had 
we started with the minimal identity 

1 + thKijOiOj = 

mj=0,l rriij=0,l 

we would have had to specify an ordering and to keep a detailed bookkeeping of 
the minus signs resulting from the rearrangement of the ipj. In this process, we 
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obtain non-local factors of the type (— l) n »J nfcl where all the indices k and I are 
distinct (or similar expressions with rriij or m^). This prevents us from performing 
the sums locally as in the above proof. However, we can use the decomposition 



2 



Each of the terms of the r.h.s can be decomposed into local factors. Each time this 
decomposition is used, the number of terms in the partition function is multiplied 
by four. We found that the number of terms in the partition function is at least 

(Af-3)(iV-4) rgi 

2 e . This result is in agreement with a remark of Kasteleyn and its proof 
in the present formulation will only be outlined. 

The lower bound on the number of terms can be obtained by establishing a 
correspondence between the non-local factors and the edge-crossings of a planar 
projection (depending on the ordering) of the complete graph' 51 with N vertices. We 
consider a surface on which the complete graph can be embedded without edge- 
crossing. From Euler theorem, the genus of this surface is at least C^^ 3 )^— 4 ) ; the 
inequality being saturated if all the faces were triangles. We then note that if the 
riij and n^i in (43) were replaced by sums, the identity would still be valid. At 
least one non-local factor of this type corresponds to the projection of each handle 
of the surface considered. This concludes the outlined proof. 



7. Conclusions 



The duality transformation and the representations in terms of Grassmann 
variables have enlighted our understanding of the nearest neighbor Ising models in 
various dimensions. We have extended these reformulations and their immediate 
applications to the case where an arbitrary number of Ising spins have a strictly 
negative interaction energy associated to any pair of spin. No other assumption 
has been made on the couplings. 
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These results apply to the hierarchical models. For these models, the renormal- 
ization group transformation is very simple and can be handled satisfactorily with 
several methods 121 when the order variables are used. The study of the continuum 
limit of the reformulations presented here using the renormalization group method 
is a challenging problem' 101 . This study is expected to shed a new light on questions 
related to the 3-dimensional nearest neighbors Ising model where the hierarchical 
approximation is rather good, as far as the critical behavior is concerned. In par- 
ticular, we would like to understand the continuum limit of the reformulation as 
"sums over equipped surfaces" proposed by Dotsenko and Polyakov. 
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and C. Zachos for valuable discussions. This paper is dedicated to the memory of 
our friends and colleagues Dwight Nicholson, Bob Smith, Chris Goertz and Linhua 
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